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A model based on frustrated total reflection of waves in a multilayered 
medium has been developed to analyze the crosstalk between multimode 
optical fibers. Kappa (k), the parameter indicating the power distribution 
among the modes of the fiber, and fiber cladding thickness play very im- 
portant roles in determining the crosstalk isolation between fibers. Sig- 
nificant but less dominant effects on crosstalk are due to variations in fiber 
numerical aperture, length, and transmitting wavelength. 

I. INTRODUCTION 

Crosstalk between communications circuits has long been a problem 
concerning engineers and designers in the telecommunications in- 
dustry. In a digital communications system using optical fibers as a 
transmission medium, at least a 30-dB signal-to-crosstalk ratio will 
be needed. In those applications where a lossy jacket around the fiber 
is undesirable, proper design of optical fiber cables and circuits requires 
an understanding of the parameters that control crosstalk between the 
fibers. 

A number of models can be found in the literature to describe cross- 
talk between optical fibers. 1-8 Maxwell's equations are the usual start- 
ing point for these models, and a field theory approach describing the 
coupling coefficients between individual modes results. Each of these 
models provides insight into a possible mechanism for describing cross- 
talk between optical fibers, and each is operationally useful for cal- 
culating crosstalk between single-mode fibers or guides with a small 
number of propagating modes. A very interesting model describing 
crosstalk resulting from scattering from a rough core-cladding inter- 
face is also described in the literature. 9 

In this paper, a meridional quasi-ray tracing procedure is used to 
describe crosstalk between highly multimoded optical fibers. This 
approach is an extension of work developed by H. P. Yuen 10 and N. S. 
Kapany. 11-14 The mechanism for crosstalk coupling of energy between 
fibers is that of frustrated total internal reflection of waves in a multi- 
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layered medium. Integral expressions for crosstalk and transmitted 
power are developed in terms of the geometry of the system, the ab- 
sorption loss of the fiber cores, the transmission and reflection coeffi- 
cients at the core-cladding interface, and the energy distribution at the 
launching end of the fiber. It is assumed in this work that the propagat- 
ing modes within a fiber are uncoupled. The error resulting from this 
assumption should be small when the excitation is chosen to coincide 
with the equilibrium energy distribution and the coupling between the 
fibers is weak. A computer evaluation has been made of the integral 
expressions for crosstalk and transmitted power. Results of a study 
showing the functional relationship between far-end equal-level cross- 
talk (fext) and cladding thickness, mode energy distribution, length, 
numerical aperture, and wavelength are included in this paper, along 
with a discussion of future work in this area. 

II. DERIVATION OF GENERAL TRANSMISSION AND CROSSTALK FORMULAS 

For the system of fibers shown in Fig. 1, we present a general deriva- 
tion of the transmission of energy within a fiber and the crosstalk be- 
havior between fibers. We assume that each fiber in an assembly is 
excited by a source that focuses its power on the center of the entrance 
end of a fiber, exciting meridional rays as shown in Fig. 2. The ray 
incident at an angle o to the axis of the fiber is refracted into the fiber 
at an angle which can be simply related to O by the law of refraction. 
Consider a distribution of input radiant intensity or input power per 
unit solid angle F (d, <f>) for some polar angle defined on the entrance 
surface. A certain fraction of the power, T (d, <f>), will be transmitted 
into the fiber because of refraction at the entrance surface. Within the 
incremental solid angle, sin dddd<j>, the power coupled into the fiber, 
dP(6, <f>), can be written as 

dP(0, <f>) = F(d, 0) sin 6d6d<t> = F (6, <}>)T (e, <f>) sin eded<f> (1) 

or 

(<?&!>) ., (M >. (1 a) 

\ (III / entrance 

We assume here that an absorption coefficient a per unit length is 
defined for the fibers. The coefficient a takes into account both bulk 
absorption loss and scattering loss. 

The power flux just prior to hitting the core-cladding interface of the 
fiber for the first time is given by 

{%) Pl =t ~ a{d/2) ^ 6F{d '^' (2) 

where d is the fiber core diameter. 
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Fig. 1 — System of optical fibers. 



The power reflected from the first reflection at the core-cladding 
interface can then be written as 



(^) fli =6-"< d/2)CflC ^(0 ) 0)F(9,<^). 



(3) 



The power transmitted into a neighboring fiber from the first reflection 
at the interface is given by : 



jjjH - e-°w™ B T(e,<{>)F(e,<t>). 



(4) 



OPTICAL 
SOURCE 




Fig. 2 — Meridional ray fiber excitation. 
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Fig. 3 — Ray tracing of transmitted and crosstalk rays in adjacent optical fibers. 

The next section shows how R(d,tj>), the reflection coefficient, and 
T{B, <f>), the transmission coefficient, are calculated. T(d, <f>) in this 
study is, in effect, the mechanism for obtaining crosstalk between fibers. 
Crosstalk is shown to be caused by frustrated total reflection of plane 
waves at the interface of a multilayered medium. The reflected rays 
continue their propagation down the fiber until they hit the wall again. 
A fraction of energy is then tunneled as crosstalk to the neighboring 
fibers, and another fraction is reflected. This process is repeated until 
the rays reach the exit end of the fiber, as illustrated in Fig. 3. Let us 
make the following definitions : 

(-j— ) = incident power density at the nth. hit 
afl / p n 

(-=— J = transmitted power density at the nth hit 
aSi Jt„ 

(jp- ) = reflected power density at the nth hit. 
du / R n 

The following power expressions for the nth hit clearly hold. 

(£),-(£),™ 

(dP\ ^(dP\ 
\ dfi ) Pn+l \ dQ, ) Bn 



ad esc 8 



(5) 
(6) 

(7) 
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For a ray of angle relative to the fiber axis, the total number of 
hits M is the largest integer smaller than (L/d)(tan 0) + \. 



M = 



^tan0+fl (8) 



and the total path length is L sec 0. The output power distribution at 
the exit end is therefore 

(tt) =e- aLB ° ee R M (e,<j>)F{6,<}>), (9) 

\ au /output 
so that the total output power of the fiber is 

r2r /-A max 

Pout P ut= / / smdF(d,<j>)R M (e,<f>)e- aLBece T (6,<l>)ddd<t,. (10) 

Jo Jo 

If the exit end of the fiber is not matched to the surrounding medium, 
there is a transmission factor, T (d, <j>) ^ 1. From eqs. (2), (5), (6), 
and (7), we obtain 

( % ) T = e ~ ad( ' l ~ i)CSC W> <t>)l n - l T(d, <j>)F(d, 0). (11) 

Let each crosstalk ray at the nth. hit suffer a further attenuation 
A' (6, 4>, L, n) before it reaches the exit end. The total crosstalk output 
power is then given by 



P xt = J'' /"""sin OF \6, 4>)T(d, 0) 

M 



X E lR(d, 0)]"- 1 ^'(0, 4>, L, n)e- adln - i)cece ddd<t>. (12) 

n=l 

Note that Mis a function of 0. Equations (10) and (12) provide the 
general formula for signal and crosstalk output power at the end of the 
fiber as a function of the system parameters. For a specific model of 
the crosstalk transfer, it is necessary to provide explicit expressions for 
F(e, <f>), R(B, <t>), T(6, 0), and A'(d, <j>, L, n). A very simple model is 
discussed later in this paper. 

That quantity of particular interest in system design is the signal- 
to-crosstalk power ratio defined by 



\NjxT 



- -^tput. (13) 

l XT 

Far-end equal-level crosstalk in decibels is then defined as 

fext = 10 log 10 ^»- (14) 

"XT 
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III. FRUSTRATED TOTAL INTERNAL REFLECTION IN A MULTILAYERED 
MEDIUM— A CROSSTALK MECHANISM 

In this section, a discussion is given of frustrated total reflection of 
plane waves in a multilayer medium. This mechanism is the one re- 
sponsible for crosstalk in our model and serves as the basis for our 
crosstalk power transmission calculation. 

Consider the usual representation of total internal reflection shown 
in Fig. 4. The ray is completely reflected when ip exceeds the critical 
angle, yp c = sin -1 (n 2 /wi). The field amplitude decays exponentially 
in the optically rarer medium. This picture is quite different when we 
look at the energy flow of the waves more closely. 16-17 The Poynting 
vector in the rarer medium is by no means zero ; only its time average 
vanishes. This energy flow is depicted in Fig. 5, which also shows that 
a displacement A exists, because of the Goos-Hanchen shift, 15,18 be- 
tween the incident and reflected beam. For plane waves, this shift is 
readily calculated. 16 It increases with decreasing angle of incidence $, 
as expected. The Poynting vector also makes a deeper penetration 
when \f/ is smaller. 

With this picture we can see that, if a third medium of refractive 
index, n\, is brought into close proximity to the first surface from below, 
as shown in Fig. 6, some energy will be trapped by it as the incident 
rays make their penetration. This trapped energy will then propagate 




Fig. 4 — Total internal reflection from plane dielectric boundary. 
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#j > ^2 > ^c 



E A 3 

Fig. 5 — Displacement of ray — the Goos-Hanchen shift. 

into the third medium. This phenomenon of energy transmission 
through a lower refractive index slab for incidence greater than the 
critical angle is called frustrated total internal reflection. It is a wave 
phenomenon outside the domain of pure geometrical optics and is 
exactly analogous to the tunneling of quantum mechanical particles 
through a potential barrier, which is classically forbidden. 

Consider now three homogeneous media that might represent the 
core, claddings, and core of adjacent fibers and model them by two half 
spaces of refractive index, n h separated by a third medium, n 2 , of 
thickness, t. We have three homogeneous media for which the trans- 
mission coefficient can be obtained by matching boundary conditions 
directly. The result, for plane wave incidence, is well known in the 




Fig. 6 — Frustrated total internal reflection in a multilayered dielectric medium. 
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literature 16 and has been derived in the appendix for the N media case 
and specialized to the three-media case discussed here. We assume that, 
in the vicinity of the core-cladding interface, the signal can be repre- 
sented as a plane wave and quote the result for the lossless case for the 
geometry and notation shown in Fig. 7. 
Let 



*-[(=)'-*- \T 



(15) 



The transmission coefficient, T(B,t), which can be defined as the 
ratio of the transmitted power flux to the incident power flux, is 
given by 

r(*,*) = [ili 

where 

Ai = the fraction of incident power polarized perpendicular to the 

plane of incidence. 
A 2 = the fraction of incident power polarized parallel to the plane 

of incidence. 

(17) 



cosh 2 + [(n 2 7 2 - m cos 2 0i)/2nin 2 cos 0ry] 2 sinh 2 J 
cosh 2 + [(n| cos 2 B x - y2n x )/2 cos dimn^J sin 2 h& J ' 



2ir . 

" T n2ty - 



Formula (16) is written in terms of d , the incident angle at the input 
of the fiber, and is used in the calculation of crosstalk power later in this 
paper. 





ni > n 2 



s 



s 




Fig. 7 — Geometry for calculation of transmission and reflection coefficients in a 
single-layered three-medium dielectric. 
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IV. CROSSTALK FORMULA ADAPTED TO CIRCULAR GEOMETRY 

We assume now that the only crosstalk that is modeled here is that 
which is due to frustrated total reflection. The geometry of the rays 
near the fiber wall at a hit is illustrated in Fig. 8, which is a cross-sec- 
tional view of a fiber and one of its neighbors. Here, <t> is the polar 
angle which also serves to define the length t. Using simple geometry 
we see that 

t = c cos <f> - d/2 - [(d/2) 2 - c 2 sin 2 «*>]* (18) 

and that the angle <t> m for which the line ON is tangent to the circle 
O'is 

m = sin- 1 d/2c. (19) 

Let E denote the center point of the entrance end of the fiber. Then a 
meridional ray ES hitting the wall of the fiber lies in the plane of inci- 




Fig. 8 — Geometry of two adjacent round optical fibers. 
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dence EOS, which cuts the fiber 0' at a line passing through S' parallel 
to the axis of the fiber. The transmission coefficient is calculated in 
the form of plane wave tunneling through a slab of thickness t, as in 
Fig. 7. When the medium between is homogeneous and lossless, eq. (16) 
can be used. 

The transmitted ray will propagate as a skew ray and emerge from 
a point a number of wavelengths from S' depending on 0, as is clear 
from Figs. 5 and 7. Since the displacement A is small compared to the 
path length between bounces and very small compared to the length 
of fiber we are interested in, we can essentially regard the ray to emerge 
from S' for the purpose of calculating the path length of the trans- 
mitted ray. 

Skew ray analysis shows 12 that the length of the skew ray crosstalk 
path lengths are also given by L sec 6. We can now calculate the ex- 
pressions for Pxt and P ut P ut- From the above discussion, substitution 
into the crosstalk integral, eq. (12), yields 

A' (n) e~ ad{n ~ i) cace = e -aL8ec9 (20) 

and 

M M 

E LR(B, 0)] n - 1 A , (n)e- ad(n - i)osoe = e^ 8609 £ R(0, 0)"- 1 

-W^*f} (21) 

The crosstalk integral simplifies to 

Pxt = f m " def 2 ' d<f> sin OF (6, tfir^'Tid, 0) 

Making a change of variable to t with 



d<}> c 2 - (d/2)* - (t + rf/2) 2 

dt " (t + d/2){(c 2 - t*)i(t + d) 2 - c 2 ]}* ' 



(23) 



we can then write 

Pxt = N j Bl " dej^ dt(^\ sin 6F{6, fle-*""^*, *) 

x[L=^f ], M 

where N = 2 times the number of fibers adjacent to the excited fiber. 
26 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1975 



V. SIMPLIFIED CROSSTALK FORMULAS USED FOR INITIAL 
COMPUTER STUDY 

A specific model requires a knowledge of F(6, <j>), R(d, <f>), and T(6, <f>) 
in eqs. (24) and (10) to calculate the crosstalk between fibers. As a 
first-order approximation, the following assumptions were made : 

(i) The interfiber medium is lossless and homogeneous so that eq. 
(16) is valid as well as: 

(it) T(d, 0) + R(6, *) = 1. (25) 

(in) The input angular power distribution of the fiber is a gaussian 
function of the form 

F(0, 0) = e-^', (26) 

where k is a parameter that is a measure of the width of the 
beam and also an indication of how the power is distributed 
between the modes of the fiber. 6 C is the critical angle of the 
fiber. 
(iv) The air spaces between the fibers are replaced by cladding 
material for the purposes of calculating the transmission co- 
efficient, T(d,4>). A partial check of this assumption showed 
that the primary contribution to crosstalk occurred for small 
values of <f>. Replacing the air spaces by cladding material 
caused a maximum error in crosstalk of 4.5 dB. 

Under the assumptions mentioned above, the crosstalk integral, 
(24), becomes 

X sin 9iT#*»*-*' m "{l - [1 - T(d, t)] M }. (27) 
The total output power of the transmitting fiber becomes 

Output = J'' d<f>[ t "*emOe- 9/ "* t e- aL '" 9 tl ~ T{e,<t>)~\ M dB. (28) 

VI. COMPUTER STUDY— SUMMARY OF RESULTS 

A computer program was written and the integrals (27) and (28) 
were evaluated for typical fiber parameters. A number of studies were 
made to determine how fext [eq. (14)] varies as a function of cladding 
thickness, numerical aperture, length, k (kappa), and wavelength. 
Table I is a guide that defines the parameters and relates the variables 
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Table I 



Far-End Crosstalk 
as a Function of 


Figure 
Number 


Parameters Held Constant 


Range of 

Independent 

Variable 


d Core Diam 
c Cladding Diam 


9 


NA = 0.10, d = 25.4 urn, L = I 

km, X = 0.6328 /xm 
a = 20 dB/km, k = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.3 to 0.7 
c 


d 
c 


10 


NA = 0.15, d = 25.4 /*m, L = 1 

km, X = 0.6328 /xm 
a = 20 dB/km, k = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.3 to 0.7 
c 


d 
c 


11 


NA = 0.20, d = 25.4 M m, L = 1 

km, X = 0.6328 Mm 
a = 20 dB/km, k = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.3 to 0.7 
c 


d 

c 


12 


NA = 0.10, d = 50.8 Mm, L = 1 

km, X = 0.6328 Mm 
o = 20 dB/km, « = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.5 to 0.9 
c 


d 

c 


13 


NA = 0.15, d = 50.8 Mm, L = 1 

km, X = 0.6328 Mm 
a = 20 dB/km, k = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.5 to 0.9 

c 


d 
c 


14 


NA = 0.20, d = 50.8 Mm, L = 1 

km, X = 0.6328 Mm 
a = 20 dB/km, k - 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.5 to 0.9 
c 


d 

c 


15 


NA = 0.10, d = 76.2 fim, L = 1 

km, X = 0.6328 Mm 
a = 20 dB/km, k = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.5 to 0.9 
c 


d 

c 


16 


NA = 0.15, d = 76.2 Mm, L - 1 

km, X = 0.6328 Mm 
a = 20 dB/km, « = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.5 to 0.9 

c 


d 

c 


17 


NA = 0.20, d = 76.2 nm, L = 1 

km, X = 0.6328 Mm 
a = 20 dB/km, k = 0.25, 0.35, 

0.4, 0.5, 10.0 


- - 0.5 to 0.9 

c 


Numerical aperture 


18 


d = 50.8 Mm, c = 85 Mm, a = 20 
dB/km, L = 1 km, X = 0.6328 
Mm, k = 0.25, 0.35, 0.40, 0.50, 
10.0 


NA = 0.05 to 
0.30 


Fiber length 


19 


d = 50.8 Mm, c = 85 Mm, a = 20 

dB/km, NA = 0.15 
X = 0.6328 Mm, « = 0.25, 0.35, 

0.4, 0.5, 10.0 


L = 100 m to 
5 km 


k (Kappa) 


20 


L = 1 km, d = 50.8 Mm, c = 85 

Mm, a = 20 dB/km 
NA = 0.15, X = 0.6328 Mm 


k = 0.1 to 1.0 


Wavelength 


21 


L = 1 km, d = 50.8 Mm, c = 85 

Mm, a = 20 dB/km 
NA = 0.15, k = 0.25, 0.35, 0.4, 

0.5, 10.0 


X = 0.6328 to 
1.06 Mm 
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09 



Fig. 9— fbxt vs. d/C (cladding thickness), d = 25.4 ixm, NA = 0.10, L = 1 km, 
X = 0.6328 /im, a = 20 dB/kra. 



in the study to the figure numbers appearing in this paper. Figures 9 
to 11 show, for numerical apertures of 0.10, 0.15, and 0.20, respectively, 
the relationships between fext and cladding thickness for a 1-mil fiber 
core diameter. For the same numerical apertures, Figs. 12 to 14 and 
15 to 17 show this relationship for 2- and 3-mil core diameters. 
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Fig. 10— pext vs. d/C (cladding thickness), d = 25 A nm, NA = 0.15, L - 1 km, 
X = 0.6328 Mm, a = 20 dB/km. 
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Fig. 11— fext vs. d/C (cladding thickness), d = 25.4 nm, NA = 0.20, L = 1 km 
X = 0.6328 urn, a = 20 dB/km. ' 
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Fig. 12— pext vs. d/C (cladding thickness), d = 50.8 A»m, NA = 0.10, L = 1 km, 
X = 0.6328 fim, a = 20 dB/km. 
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Fig. 13 — fbxt vs. d/C (cladding thickness), d = 50.8 urn, NA = 0.15, L = 1 km, 
X = 0.6328 Mm, a = 20 dB/km. 



OPTICAL FIBER CROSSTALK 



33 



100 
90 
80 

70 1— 



40 



30 



10 



0.3 




L=1 km 
X=0.6328/xm 
a=20dB/km 



0.5 0.6 0.7 0.8 0.9 

d/C 
(25.4) (16.9) (10.9) (6.35) (2.82) 

CLADDING THICKNESS IN MICROMETERS 



Fig. 14— fext vs. d/C (cladding thickness), d = 50.8 Mm, NA = 0.20, L = 1 km, 
X = 0.6328 Mm, a = 20 dB/km. 
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Fig. 15— fext vs. d/C (cladding thickness), d = 76.2 urn, NA = 0.10, L = 1 km, 
A = 0.6328 y.m, a = 20 dB/km. 
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Fig. 16— fext vs. d/C (cladding thickness), d = 76.2 urn, NA = 0.15, L = 1 km, 
X = 0.6328 fxm, a = 20 dB/km. 
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Fig. 17— fext vs. d/C (cladding thickness), d = 76.2 Mm, NA = 0.20, L = 1 km, 
X = 0.6328 urn, a = 20 dB/km. 
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For practical claddings greater than 12 microns in thickness, in- 
creasing the cladding thickness will improve crosstalk isolation by 
approximately 0.8 dB/micron. Cladding thickness is an important 
design parameter for improving crosstalk isolation between optical 
fibers. If we attempted to eliminate crosstalk by coating the cladding 
with an opaque substance, cladding thickness would also play an 
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Fig. 18 — fext vs. numerical aperture, d = 50.8 Mm, C = 85 Mm, a = 20 dB/km, 
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important role in determining the amount of loss suffered by the trans- 
mitted energy in the core owing to the lossy coating. A cross-check of 
Figs. 9 to 17 shows a very weak dependence of fext on core diameter. 
A fiber core diameter of 50.8 nm and cladding diameter of 85 /xm was 
chosen to represent a typical fiber geometry in the remainder of this 
study. Figure 18 shows the relationship between fext and numerical 
aperture with a fixed cladding refractive index of 1.458. For a given 
kappa, crosstalk isolation improves by approximately 6 dB by doubling 
the numerical aperture in the range from 0.10 to 0.30. For an NA 
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= 0.15, crosstalk was evaluated as a function of length as shown in 
Fig. 19. For a fixed kappa, crosstalk isolation decreased by approxi- 
mately 2 dB per decade of length. The effect of kappa on crosstalk is 
shown in Fig. 20. Kappa is an indication of the distribution of power 
among the modes of the fiber. As kappa decreases in value, the energy 
distributed in the lower-order modes of a fiber increases. A knowledge 
of the steady-state mode distribution should enable us to estimate an 
effective k for use in this model. As illustrated in Fig. 20, k is an ex- 
tremely important parameter in determining the crosstalk between 
fibers. For kappas less than 0.5, crosstalk isolation is greater than 40 
dB for a kilometer length of fiber. For kappas less than 0.40, crosstalk 
isolation is greater than 50 dB. The final study, illustrated in Fig. 21, 
shows the relationship between fext and wavelength. For a given 
kappa, doubling the wavelength of the transmitting signal decreases 
crosstalk isolation by approximately 7 dB. It is envisioned that typical 
fiber losses in a future optical transmission system will be less than 
the 20-dB/km loss reported here. We will extend, in the future, our 
computer study to include fiber losses ranging from 2 to 20 dB/km, 
but do not expect this loss parameter to change the general conclusions 
drawn in this paper. 

To determine if the primary mechanism for crosstalk in optical 
fibers is frustrated total reflection of waves in a multilayered medium, 
a crosstalk experiment on a long length of parallel fibers must be per- 
formed. In this experiment, crosstalk and kappa should be measured 
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as a function of length at a number of different wavelengths. This type 
of experiment will enable us to check many aspects of the model. 

APPENDIX A 

Calculation of Transmission Coefficients in a Multilayered Medium 

The general expression for the transmission coefficients in a multi- 
layered dielectric medium is well known in the literature 151718 and is 
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presented here. This expression is simplified for the three-medium case 
and, ultimately, eq. (16) in the text is derived. 

Consider the geometry shown in Fig. 22. Let us suppose that between 
two semi-infinite media, denoted by 1 and n + 1, there are n — 1 
layers of dielectric material denoted by 2, 3, • • • , n. Let a plane wave 
be incident on the last layer at an angle of incidence 0„+i and let the 
plane of incidence be the x - z plane. As a result of multiple reflec- 
tions at the boundaries of the layers, two waves exist in each medium 
with the exception of medium 1. Our problem will be to determine the 
amplitude of the transmitted wave in medium n+1 and hence the 
transmission coefficient. The following notation will be used. 

Zj = the coordinate of the boundary between the jth and ( j -f l)st 

layers. 
dj = Zj — Zf-i = the thickness of the jth layer. 
kj = (2ir/X)(ny) = the wave number in the jth medium. 
a, = kj cos 9, = z component of the wave vector in the jth layer. 
<f>j = otjdj = the phase change in the jth medium. 
Zj = the self-impedance of the jth layer. 

Z\ n = the input impedance looking into the jth medium from the 
j + 1 medium. 

The electric and magnetic fields in the jth medium can then be 
written as 



E jv = Ajexp l-ictjiZ - Zj- X )~\ + £ y exp [wtj(Z - Z/_i)] 



(29) 



H* = ± [Ajex V l-iaj(Z - Zj^n ~ Bjexp \iutf - Z^)]}. (30) 
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Fig. 22 — Geometry used for calculation of transmission coefficients in multilayered 
dielectric media. 
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The x and t dependency in this case is omitted for the sake of brevity, 
but assumes the general form : 

exp i(k n+ ix sin n+ i — tot). 

Aj and B h B x = are the amplitudes of the incident and reflected 
waves in the jth medium. The amplitude -A„+i of the incident wave is 
assumed to be known. To obtain the transmission coefficient of interest, 

_ Ai _ AjA 2 A„-i A„ . > 

Tn+1 A n+ i A 2 A S An A n +i 

We can write 2n boundary equations for the tangential components of 
the fields and solve these equations for A h A 2 , ■••, An+i', B\ t B 2 , • • •, 
B n+ i. When these coefficients are known, the transmission and reflec- 
tion coefficients for the multilayered medium are obtained. 

A second approach, and one that develops an iterative scheme more 
suitable to a digital computer, describes the transmission coefficient in 
terms of a generalized input impedance. 19 This is the approach that 
will be followed here. It is straightforward to derive, and it is shown in 
the literature that : 

i k , = / iJ TT5 exp(i *' ) - (32) 

Substituting into eq. (31) yields 

rn + i = II f' + _5 exp fa), (33) 

with di = 0, where 

lij — tZiia tan <f>j 

Specializing eq. (33) for the three-medium case discussed in the text, 
we can obtain, with some algebraic manipulation, the following 
formula : 

, = *Mi ,»*> 

72 (Z, - Z 2 )(Z 2 - Z 3 )e**- + (Zi + Z 2 ) (Z, + Z 3 )r^' K60) 

For electric fields parallel and perpendicular to the plane of incidence, 
eq. (35) will take different forms. For E x , 

Zi = Z ° Q = Z, (36) 

n\ cos di 

Z 2 = —^— = A , (37) 

n 2 cos 6 2 m 2 y 
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where 

*-[(2) i -*- i r (38) 

Z = ^ (39) 

/3 = k 2 d 2 y. (40) 

After some algebraic manipulation, 

_ Kj cosh fl - iKiKi sinh <3 ^ 
T2i ~ i^ cosh 2 + #l sinh 2 ' 

Ki = 2ni7i 2 cos 0i7 (42) 

K 2 = n! 7 2 - n\ cos 2 6 h (43) 



where 



for^i 



= Z t COTfli = Za (44) 

= Z O COS0 2 = tZoT (45) 

72.2 n 2 



and 2Tsi, becomes 

K\ cosh fl + tgjgs sinh fl 
T2 " " K 2 cosh 2 + K§ sinh 2 ]8 



(46) 



where 

K a = wl cos 2 0i - 7 2 w?- (47) 

Utilizing eqs. (41) and (46), we can obtain the power transmission co- 
efficients directly : 

= cosh 2 18 + [(nl7 2 - nf cos 2 di)/2iiin 2 cos 0ry] 2 sinh 2 

^2„ = |T 2lI | 2 

1 . (49) 

cosh 2 + [(nl cos 2 0i - y 2 n\) /2 cos 0in 1 n 2 7] 2 sinh 2 

Equations (48) and (49) were used to obtain eq. (17) in the text of this 
paper. 
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